Ground State Properties of Spin-Orbit Coupled Bose Gases for Arbitrary Interactions 



(N 

o 



Renyuan Liao 1 ' 2 and Wu-Ming Liu 1 
1 National Laboratory for Condensed Matter Physics, Institute of Physics, 
Chinese Academy of Sciences, Beijing 100190, China and 
2 College of Physics and Energy, Fujian Normal University, Fuzhou 350108, China 

(Dated: September 14, 2012) 

We study spin-orbit coupled (SOC) Bose gases with arbitrary interspecies interaction. Besides 
at a critical interaction our results carry over to the recent publication [PRL 109,025301 (2012)], 
we identify various new features arising from the interplay of SOC and interspecies interaction, 
including a roton minimum in the excitation spectrum and dual effects of SOC on ground state 
energies depending on interspecies interactions. Counterintuitively, we find that at low interspecies 
interaction the SOC stabilizes the system by suppressing the quantum depletion. We show that 
the static structure factor is robust with respect to the SOC in the phase space where time-reversal 
symmetry is preserved. 
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The pioneering experimental realization of synthetic 
gauge field and spin-orbit coupling (SOC) [lj-Q provides 
fascinating opportunities to explore quantum many-body 
systems of ultracold atomic gases. Synthetic gauge po- 
tential is defining a new dimension for simulating real 
materials with cold atoms. The engineered spin-orbit 
coupling in a neutral atomic Bose-Einstein condensate 
was accomplished by dressing two atomic spin states with 
a pair of lasers. In electronic systems, spin-orbit coupling 
is crucial for quantum spin Hall effects and topological in- 
sulators Q , which has captured a great deal of attention 
in condensed matter community. In particular, spin-orbit 
coupled bosons does not have an analogy in conventional 
condensed matter systems, resulting in the emergence of 
many novel quantum phases such as striped supcrfluid 
phase 0-Q and half vortex phase l^- 14 1 . 

Inspired by experimental achievements, theorists have 
been paying vast attention to SOC Bose gases in many 
aspects, including the fluctuations effects 
quasiparticle excitation spectrum 18 , Hij . 
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These theoretical works focus mainly 
on two types of spin-orbit coupling: one is the SOC 
with equal Rashba and Dresselhaus couplings firstly re- 
alized at NIST l], the other is the Rashba SOC which is 
theoretically more complicated and_ appealing. For the 



coupling 22|. In addition, we find that by tuning the 
interspecies interaction to a critical one, a second gap- 
less mode emerges, consistent with Ozawa and Baym's 
result 20]. Through inspecting the infrared behaviors 
around the gapless modes, we are able to conclude that 
the Bose-Einstein-Condensates (BEC) is stable, in spite 
of the fact that in the absence of interactions there is 
no BEC for SOC Bose gases in three dimensions due to 
constant density of states at low energy. Secondly, by 
studying the ground state energy, we find that at low 
interspecies coupling, the ground state energy actually 
increases with SOC, in contrast to that of large inter- 
species coupling. Thirdly, we clarify the effects of SOC 
on the quantum depletion. At low interspecies coupling, 
we show that the SOC actually suppresses the quantum 
depiction, differing from a reported result 2l|. Finally, 



we find that the experimentally important quantity, the 
static structure factor is robust with respect to SOC in 
the parameter space where the time-reversal symmetry 
is preserved. 

We consider three-dimensional homogeneous two- 
component Bose gases with an isotropic in-plane Rashba 
spin-orbit coupling, described by the following grand 
canonical Hamiltonian: 
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Rashba-type SOC, these works 

dedicated to the critical case where the intraspecies and 
interspecies interactions are the same, leaving the physics 
of the interplay between spin-orbit coupling and inter- 
species interaction, which are both experimentally rele- 
vant and theoretically interesting, largely intact. 

In this work, by monitoring the effects of SOC and 
interspecies interaction, we have identified several new 
features and clarified several important issues regard- 
ing SOC Bose gases. Firstly, by examining the excita- 
tion spectrum, we find a roton minimum, a result also 
found for the SOC with equal Rashba and Dresselhaus 
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where <?i is the intraspecies interaction, g\ 2 is the inter- 
species interaction, and the spin-orbit coupling is R — 
\{P X — iPy), with A being the coupling strength. 

Within the framework of imaginary-time field integral, 
the partition function of the system in the momentum 
space can be cast as Z = J ^> q(7 ] exp (— S[ipi a , tpqo-]) 
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with the action 
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Here, we have defined the free particle energy disper- 
sion £ q = q 2 — fi (we have set H = 2m = ks = 1) 
and i? q = \q±e~ tVc >, with ^ being the magnitude of 
in-plane momentum, and tpq — Arg(q Xl q y ). The chem- 
ical potential /i is introduced to fix the total number of 
the particles, and V is the volume of the system. For 
a non-interacting system, the Hamiltonian is diagonal- 
ized in helicity basis with a dispersion — £ q ± Xq±. 
The lowest energy states are infinitely degenerate, sit- 
ting at the circular ring in the plane q z = confined 

by <lx + Qy — (§) m the momentum space. Mean-field 
study find that the plane wave (PW) phase exists 
for <?i > <?i2 and the striped phase exists for 171 < 312. 



The PW phase is characterized by a condensation at 
a single momentum state, while the striped phase is a 
coherent superposition of two condensates at two op- 
posite momenta. Here we restrict ourself to studying 
the PW phase, and further assume that the condensa- 
tion occurs at momentum k = (A/2,0,0). We sepa- 
rate a Bose field into a mean-field part and a fluctuating 
part: Vqo- = 4>ca5qR + </>qer- Without loss of general- 
ity, the condensate wavefunction in spin space could be 



chosen as $(x) 



y = vT(i,-i) 
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that \<t>ca\ = \fn^a = v no/2, with no being the total 
condensate density. Retaining terms of the zeroth and 
quadratic orders in the fluctuating fields, we rewrite the 
effective action as S e f / = S c + S g . Here the mean-field 
action is S c = /3VJ2 C [(£ K - Ari)n C(T + (g x + g\ 2 )n 2 ca \ . 
Saddle point conditions 5S c /8n ca = leads to fj, = 

~X + + 9u)no and S c = ~^f(gi + g^hl- 
By defining a four-dimensional column vector $ q = 

(<rWqt)0;?+q-;>0!-qt>$£-q4.)' W6 CEn bring th<3 S auSsian 

action into a compact form S g — J2 q ^qp~ 1 ^q — 
/3^ q e_ q , where q = (q, iw n ) with w n — 2i:n/j3 be- 
ing the bosonic Matsubara frequencies, and the inverse 
Green's function Q^ 1 (q,iw n ) defined as 
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where e q = ^ +q + (2g x + g 12 )n = %- + q 2 + Xq x + gin Q 
and i? q = X(k + q)x. Throughout our calculation, we 
will choose g\ n as a basic energy scale and yfg\n^ as 
the corresponding momentum scale. We define a dimen- 
sionless parameter -q = 312/31 {v € [0, 1]) to characterize 
the strength of interspecies interaction. 

The excitation spectrum provides useful insights into 
a system. It could be obtained by examining the poles 
of the Green's function. To achieve this, one proceeds by 
evaluating the determinant of <7 -1 (q, iw n ), 

DetiQ- 1 } = {iw 2 n - w 2 ) [(iw n - 2\q x f - w 2 ] - 2A 2 g 2 F, 
F = iw n (iw n - 2\q x ) + q 2 (q 2 + 2gxno + A 2 ) + 

X 2 q 2 

(gi + g 12 )n [(gi - gi 2 )n + A 2 ] ^, (4) 

w here uj 10 = y 1 q 2 [q 2 + 2(ffi + g 1 2 )n ] and uj 20 = 
{q 2 + A 2 ) [q 2 + A 2 + 2(31 - 512)710]- By requiring 
Det[g^ 1 (q, iw n )] = 0, one can obtain the excitation 
spectrum of the system. It is easily to verify that 
Det [t/ _1 (— q, — iw n )] = Det [^/ _1 (q, iw n )] , so that it 
has two excitation branches. 

In the q x -q z plane (namely q y = 0), it is straight- 
forward to analytically derive two branches of the ex- 



citation spectrum from Equation uj\ = u3\q and 

lj 2 = 2\q x + W20- For the first branch, there ex- 
ists a gapless mode at qi = (0,0,0). Expanding 
around this gapless point, one finds that the low-energy 
spectrum is phonon-like collective excitation wi(5q) ~ 

y/2( 5 i + g 12 )n a (Sq 2 + Sq 2 ) + 0(6q*). For the second 

branch, it is always gapped if 31 > 312. At the crit- 
ical point 31 = 312, a second gapless mode emerges 
at q 2 = (—A, 0,0). Expansion around this gapless 
point leads to a free particle-like anisotropic excitation 
w 2 (q 2 +<5q) ~ Sql + 5ql+g 1 n /(X 2 + 2g 1 n )Sq 2 . Judging 
from the behaviors of the low-energy spectrum around 
the two gapless modes, we conclude that there is no in- 
frared divergence for quantum depletion of the conden- 
sate in three dimensions at zero temperature as the inte- 
gral / q 2 dq sin 9 'd(p/w s (q)(s = 1,2) is finite. 

The behavior of these two branches of the excitation 
spectrum are summarized in Fig. [TJ As the spectrum en- 
joys the inversion symmetry along x-axis: ui±(q x , Qy,0) = 
^±(&) — Q-y, 0), we only plot the spectrum for the az- 
imuthal angle (p lying in [0,7r]. For intermediate inter- 
species interaction r\ = 0.5, we plot the two branches in 
panel (a) and (b) . The lower branch w_ (Fig. [1^,) accom- 
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FIG. 1. (Color online) The excitation spectrum for A = y/gino 
at zero temperature: (a) the lower branch ui- for r\ — 0.5, (b) 
the upper branch ui+ for 77 = 0.5, (c) the lower branch oi_ 
for 77 = 1.0, and (d) the upper branch w+ for ij = 1.0. Here 
V = ffia/ffi- 



FIG. 2. (Color online) The shift of the ground state energy 
AE G = E G (X) - E G (X = 0) [in units of V(gwo) 5/2 ] as a 
function of spin-orbit coupling strength A (in units of ^/gyno) 
for different interspecies interaction parameter rj: (a) 77=0, 
(b) rj = 0.3, (c) rj = 0.6, and (d) 77 = 1.0. Here 77 = 312/31- 



modates gapless excitation at zero momentum along all 
directions. For p> = 0, 7r/4, tt/2, and 37r/4, the quasipar- 
ticle energy for both branches increases monotonically 
with the momentum. Remarkably, for tp = 71", there ap- 
pears a roton minimum, signaling that the system has 
the tendency toward crystallization. For the critical in- 
terspecies interaction 77 = 1, we plot the two branches in 
panel (c) and (d). For the lower branch, the quasiparticle 
energy for both branches increases with the momentum 
for ip = 0, 7r/4, 7t/2. For ip — 37r/4, both branches have 
a local minimum at certain finite momentum. At ip = n, 
a second gapless excitation develops at q = (—A, 0, 0). 

The thermodynamic potential is given by 51 = 
— lnZ//3 = O c + fig, where f2 c = —V{g\ + g\2)rv^j2 and 
il g = ^TrlnC/ -1 — X) q e q- The thermodynamic poten- 
tial £1 possesses an ultraviolet divergence, an artifact of 
zero-range interactions, which can be removed either by 
replacing the bare interactions g\ and 512 with T-matrix 
or by subtracting counter terms [24j. At zero tempera- 
ture, the ground state energy becomes Eg = O + fJ.N, 
renormalized as 
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Here Emf — V{g\ + gi2)nl/2 is the mean-field energy, 
independent of the SOC strength A. The effects of spin- 
orbit coupling on the shift of the ground state energy 
AE G = E G (X) ~ E G (A = 0) is shown in Fig. H In the 
absence of interspecies interaction (Fig. [2^), the ground 
state energy increases with the SOC strength A. At 
rj = 0.3, AEq increases with SOC when SOC is small, 



and decreases with SOC when SOC is large. When 77 is 
sufficiently large (Fig. [2j:,d), AEg decreases with SOC 
monotonically. In the absence of the interspecies inter- 
action and SOC, we have verified that the ground state 
energy for either species recovers the Lee-Huang- Yang re- 
sult [25| for spinless weakly-interacting Bose gases with 
E G /V = ^(1 + y|^)\/na3, where a is the scattering 
length. Up to the gaussian level, the ground state en- 
ergy is universal in the sense that it depends solely on 
the gas parameter na 3 and not on microscopic details of 
the interaction potential. 

Being an intrinsic property of a BEC, the quantum 
depletion provides key information about the robustness 
of the superfluid state. The number of excited particles 
is given by 



= 2 ^2 Gn(q, iw„). 



(6) 



We show the quantum depletion of the condensate 77, ex in 
Fig. |3l In the absence of interspecies interaction (77 = 0, 
see Fig. |3k), n cx develops a global minimum at a crit- 
ical SOC. At sufficiently large interspecies interaction 
(77=0.5, 1), tt-cx increases monotonically with A, leading 
to an enhanced quantum depletion. The effect of in- 
terspecies interaction on 77, cx is shown in panel (b) for 
three typical spin-orbit coupling strengths A=0, 0.4 and 
0.8. For a fixed SOC strength , the number density of 
the excited particles increases monotonically with inter- 
species coupling parameter 77. At low interspecies inter- 
action, the SOC and interspecies interaction counteract 
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FIG. 3. (Color online) The number of excited particles n ex 
[in units of (<?i77o) ] as a function of (a) spin-orbit coupling 
strength A (in units of y/gino) for different interspecies inter- 
action parameter 77, and (b) interspecies coupling parameter 
77 for different spin-orbit coupling strength A. 



on the quantum depletion; while at sufficiently large in- 
terspecies coupling, they work cooperatively. Being pro- 
portional to (gino) 3 / 2 , is small when Bogoliubov con- 
dition gino <C 1 is satisfied, justifying our treatment. 

The static structure 5(q) probes density fluctuations 
of a system. It provides information on both the spec- 
trum of collective excitations, which could be investi- 
gated at low momentum transfer, and the momentum 
distribution, which characterizes the behavior of the sys- 
tem at high momentum transfer, where the response is 
dominated by single-particle effects. We can evaluate the 
static structure factor at the Bogoliubov level as follows: 



NS(q) = < SpUp* > 



f ££(-i) i+ ^ 



(-2q 2 )A(q,iw n ) 
DetiG- 1 ) : 
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where A(q,iw n ) = (iw n -2\q x ) 2 -(q 2 +\ 2 -\ 2 q 2 /q 2 )[q 2 + 
A 2 + 2(<7i — 512)^-0] • In the Qx — <Zz plane, namely q y = 0, 
we have 



S(q x ,qz 



No 
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This is exactly the Feynman relation [261 ] , which connects 
the static structure factor to the excitation spectrum for a 
Bose system with time-reversal symmetry. Therefore for 
q y =0, we find that the Feynman relation is preserved. 
For any other direction, numerical results suggest that 
the Feynman relation is violated. This is reasonable as 
the ground state breaks time-reversal symmetry except 
in the q x -q z plane. S(q) possesses the inversion sym- 
metry S(— q) = S'(q), and increases with the magnitude 
of the momentum, approaching unit at high momentum. 
The static structure factor 5*(q) in terms of in-plane mo- 
mentum (we choose q z — 0) is shown in Fig. |4] For the 
upper panel, we plot S'(q) at A = 3 (in units of y/gTno) 
for different interspecies interaction (r/=0, 0.5, 1). S(q) 




FIG. 4. (Color online) The static structure factor S(q) as a 
function of in-plane momentum q = (q x , q y ,0). Upper panel: 
at A — 3 (in units of ^/gTno) for different interspecies coupling 
strength (a) 77 = 0, (b) 77 = 0.5, and (c) 77 — 1.0. Lower panel: 
at 77 = 0.5 for different spin-orbit coupling strength A (d) 
A = 0, (e) A = 4, and (f) A = 8. The contours on the bottom 
plane is the projection of S(q), presented for aesthetic appeal. 



increases monotonically with in-plane momentum but de- 
creases with 77. The interspecies interaction 77 has similar 
effects on S*(q) along g x -axis and g^-axis. In the lower 
panel, we plot 5(q) at 77 = 0.5 for different SOC (A=0, 
4, 8). 5(q) along y-axis decreases with SOC, in contrast 
to the case along g x -axis, where 5(q) remains the same. 
This could be explained from Equation ([8} as 5(q) does 
not depend on A for q y = 0. 

Our predictions bear consequences for experimental 
observation. The anisotropic nature of the excitation 
spectrum could be probed by momentum-resolved pho- 
toemission spectroscopy 27]. Via measuring in situ den- 
sity distribution, one can determine the ground state en- 
ergy of the system [28| . The depletion of the condensate 
could be directly observed as a diffuse background in the 
time-of- flight images j29[. The Bragg spectroscopy can 
be employed to measure the static structure factor of the 
system [301 ] . 

In summary, we have studied spin-orbit coupled Bose 
gases with intraspecies and interspecies interactions. We 
have identified various new features arising from the in- 
terplay of the SOC and interspecies interactions. We 
hope that our work will add new excitement to the surg- 
ing field of spin-orbit coupled quantum gases. 
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